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ABSTRACT 

We present a self-similar solution to describe the propagation of a shock wave whose energy is 
deposited or lost at the front. Both of the propagation of the shock wave in a medium having a 
power-law density profile and the expansion of the medium to a vacuum after the shock breakout are 
described with a Lagrangian coordinate. The Chapman- Jouguet detonation is found to accelerate the 
medium most effectively. The results are compared with some numerical simulations in the literature. 
We derive the fractions of the deposited/lost energy at the shock front in some specific cases, which 
will be useful when applying this solution to actual phenomena. 

Subject headings: hydrodynamics — shock waves — supernovae: general 



1. INTRODUCTION 

The self-similar solution is a powerful method to de- 
scribe flows of a fluid involving strong shock waves. In 
a supernova (SN) explosion, a shock wave propagates in 
the stellar interior toward the surface. This process can 
be regarded as a self-similar problem as long as the den- 
sity distribution has a power-law profile with respect to 
the distance from the center or the surface. 

In fact, self-similar solutions for this problem 
have been discovered an d w idely kn o wn ( see, e.g., 
IZel'dovich fe Raizeri l2002f ). ISakurail (|1960l ) discov- 
ered a self-similar solution for the emergence of a 
strong shock wave from a medium with a power- 
law density distribution, which was first posed by 
iGandel'man fe Frank-KamenetskiH (fl956). This solution 
describes distributions of the velocity, pressure and den- 
sity of the matter with an Eulerian coordinate until the 
shock wave reaches the surfac e. After that, the solution 
uses a Lagrangian coordinate. iMatzner fe McKee] (|1999f ) 
invented another description for a whole evolution of this 
solution with a Lagrangian coordinate. 

However, these solutions do not completely reflect on 
the circumstances. In reality, thermal energy can be de- 
posited and/or lost across a shock front (referred to as 
the energy variation in this paper). For example, some 
nuclear reactions can take place in the region behind the 
shock front when the temperature thereof becomes suf- 
ficiently high. This happens in a detonation wave and 
might be realized in a type la supernova (SN la), which 
is thought to be a thermonuclear explosion of a white 
dwarf. The thermonuclear explosion mode l for S Ne la 
was originally proposed bv iHovle fc Fowler! ([i960), and 
several models for the propagation of the flame front have 
been discussed, for example, "deflagra tion model" (see 
llvanova et alJ fi974: Nomoto ct al. 1976), " delayed deto- 
nation model" (seelKhokhlovll991atlYamaoka et alll992t 
IWooslev fc Weaverlll994fl, and "pulsating delayed deto- 
nation model" |Khokhlovlll991b( l. In particular, the (pul- 
sating) delayed detonation models suggest that the det- 
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onation wave is generated somewhere in the outer layer 
of a white dwarf and propagates toward the surface. In 
addition, a carbon detonation model of a white dwarf 
with the mass of ~ 1 M Q was proposed (|Shigevama et al.1 
11992( 1 for a peculiar SN la, SN 1990N. The detonation 
wave propagates toward the surface also in this model. 
Furthermore, thermal energy generated behind the shock 
front can be transported through the front near the stel- 
lar surface where the optical depth is close to unity. In 
this region, photons diffuse out at velocities close to the 
speed of light, which is larger than the shock velocity. In 
other words, the ener gy of the shock wa ve can decrease by 
the radiative cooling. IChevalie"rl (|1976( l calculated the ef- 
fect of the cooling using a radiation-hydrodynamics code. 

This kind of phenomena in some type I SNe accel- 
erates a considerable amount of matter to sufficiently 
high energies to synthesize light ele ments, Li, Be, and 



B through spallation reaction s (see iFields et al. 120021 ; 
iNakamura fc Shigevamal l2004t iNakamura et all 12006( 1 . 
Thus the solution derived in this paper might be useful 
to investigate the chemical evolution of light elements in 

g alaxies. 

iBarenblatt fc Sivashinskiil (|1970f l studied the treat- 
ment of the energy variation at the shock front in the 
frame work of self-similar solutions and showed that it 
can be treated by modifying the jump conditions at the 
front to impose the condition that the energy deposition 
or loss per unit mass be proportional to the internal en- 
ergy. In real shock acceleration problem, this condition 
is rarely satisfied. Whether the energy is deposited by 
nuclear reactions or lost by radiation, it tends to have 
a characteristic energy per nucleon or a characteristic 
time scale. Thus the self-similarity of the associated flow 
might appear in very limited space and/or time if any. 

In the next s ection, we review the procedure of 
iBarenblattl (|1996l ) to treat the energ y variation. In § 3, 
we for mulate the problem following Ma tzner fc McKea 
(1999). The integrations of the governing equations are 
found to be reduced to an eigen value problem. The Ap- 
pendix gives a method to determine the eigen value. The 
results are presented in §4. Furthermore, the asymptotic 
behaviors of the velocity, pressure, and density of the 
ejecta are derived. The energy distributions are deduced 
from the exponents of these physical quantities (§5). We 
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TABLE 1 
EIGEN VALUES A FOR 7 = 4/3, 













7c 












a 


11/3 


3.0 


2.0 


1.5 


4/3 


1.2 


1.1 


1.01 


1.001 


1.0001 


3 

1.5 


1.0000 
0.5000 


0.9525 
0.4725 


0.7915 
0.4000 


0.6407 
0.3271 


0.5572 
0.2863 


0.4608 
0.2386 


0.3484 
0.1822 


0.1212 
0.06507 


0.03733 
0.02050 


0.01092 
0.006113 



conclude this paper in § 6. 

2. TREATMENT OF ENEGY VARIATION 

In this section, we review the procedure to treat the 
energy deposition or loss at the shock front introduced in 
iBarenblattl (|1996l ). The shock propagates in the ideal gas 
whose density is po . We assume that the gas ahead of the 
shock front is stationary and neglect its pressure. The gas 
behind the shock front has the velocity itf, the pressure 
Pi, and the density p?. Then the energy conservation 
law across the shock front propagating at the speed U is 
described as 

7 pf (v,[-U) 2 ~ 



Pi{u { -U) 

_7 — lpf 2 
Using the mass conservation law 

pt(ut - U) = -poll, 
and the momentum conservation law 
pi(u{ - U)v.{ + pi = 0, 

we obtain 

Pi , Uf 

2 



PoU* 



Pi(ui - U) 



q 



+ P{U{ = 0, 



(7 - i)w 

where we have added an energy flux term — p{(u{ 



(1) 



(2) 



(3) 



(4) 
U)q 



to the left hand side of equation (p}. Here the energy 
source term is denoted by q. This source term needs to 
be proportional to the internal energy to guarantee the 
self-similarity of the problem. Thus the source term is 
expressed as 

7c - 7 Pi / K s 

q =7 7T7 TT — ' w 

(7c - 1)(7 " 1) Pf 

by introducing the effective adiabatic exponent 7 C 
(Barcnblatt & Sivashinskii 1970). The physical meaning 
of 7c is as follows; (1) if 7 C > 7, energy is deposited at the 
front. (2) if 7 C = 7, the situation is same as the classical 
very intense explosion problem. (3) if j c < 7, energy 
is lost at the shock front. Using the effective adiabatic 
exponent 7 e , we transform equation (pf]) to 

„,2" 



pt(u { - U) 



Pi 



+ PiUi = 0. 



(6) 



.(7c-l)pf 2 
This expression has the same form as the classical very 
intense explosion problem except for replacing 7 with j e . 
Therefore, the Rankine-Hugoniot relations for a strong 
shock become 

2 



7e 

2 



1 



u. 



Pi 



7c + 1 



P0U 2 



Pi 



7c + 1 



7P0, 



(7) 



(8) 



(9) 



7c -1 

In other words, we can treat the energy variation by vary- 
ing 7c 



3. FORMULATION 

In this section, we derive equations describing the 
self-similar motion of t he stellar matter. Following 
iMatzner fc McKee] (|1999l) , we treat the shock emergence 
with a Lagrangian coordinate. 

3.1. Basic equations 

We first define the space coordinate x as the distance 
of the fluid element from the stellar surface and the time 
coordinate t as the time measured from the moment when 
the shock reaches the surface. The initial position of each 
fluid element is denoted by xq. The basic equations to 
describe the evolution of the density p, velocity u and 
pressure p of the fluid are formulated with the Lagrangian 
coordinate xq as 

dx _ po 
dx p ' 



(10) 



du 
PP~ 



dx 

m' 

1 dp 
Po dx ' 

z PiP7 J - 



(11) 
(12) 
(13) 

We assume that the stellar envelope has a power-law den- 
sity profile; 



Po[X) ~ 1 for x < 0, 



(14) 



where k\ and a are constants. If the envelope of the 
star is radiative, the power-law exponent is a = 3. On 
the other hand, a = 1.5 corresponds to the convective 
envelope. The time when each fluid element experiences 
the shock is denoted by to(< 0). By defining the variable 
m as 



Pa(y)dy, 



(15) 



the coordinates x, Xq, and to are labeled by m. Then, 
we define the similarity variable rj as 

il = TT-y ( 16 ) 
to(m) 

and the flow variable S(rj) as the position of a fluid ele- 
ment normalized with its initial position, i.e., 



S(r)) = 



x(m) 
x (m)' 



(17) 



Both of these variables rj(m) and S(rj) become unity 
when the shock hits the fluid element with m whereas 
•q becomes zero when the shock reaches the surface. The 
phase with < 77 < 1 corresponds to the evolution of 
shocked fluid elements till the shock breakout and the 
phase with 77 < to the expansion of the stellar matter 
to a vacuum. 
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1. — Profile of the flow variable S(tj) for a - 
= 11/3,2.0,4/3,1.2,1.01 

3.2. Expression of physical variables 



3.0,7 



Next we derive relations between the flow variable S(rj) 
and the velocity, pressure and density of the stellar mat- 
ter. We assume that the velocity of the shock front takes 
the form of some power of the original position xq of the 
fluid element; 

U = -k 2 [x {m)]-\ (18) 

where and A are constants. The latter will be deter- 
mined as the eigen value. Then we obtain the following 
relation by integrating equation (| 1 8|) with respect to time 
from t = to t = to 



in 



(19) 



fe 2 (l + A) 

On the other hand, a differentiation of both sides of equa- 
tion (|16p with respect to xq leads to 



drj = —77(1 



Xq 



(20) 



Using equations (fT71) . (l20p . and (flU)) , the density is ex- 
pressed as 

Po 



(21) 



5(77) -(a + i) v s'( v y 

where the prime denotes the derivative with respect to 
rj. Then the pressure is also expressed in terms of S(rj) 



as 



, P 

P=Pf — 
.Pi 



P0U 2 



the velocity 



x[S(rj)-(\ + i)r,S'(ri)]- 
u = (1 + X)S'U. 



(22) 



(23) 



As a result, solving equations (fT0|) through (fl3|) reduces 
to the integration of a second order ordinary differen- 
tial equation for the flow variable S(n). The method to 
determine the function S(rj) is described in Appendix. 




-100 -80 -60 -40 -20 

Fig. 2. — Profile of the velocity u/Uf for a = 3.0, 7 = 4/3, y c 
11/3,2.0,4/3,1.2,1.01 



Q. 
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Fig. 3. — Profile of the density p/pf for a = 3.0,7 = 4/3, 70 = 
11/3,2.0,4/3,1.2,1.01 

4. RESULTS 

4.1. Distributions of the physical variables 

The distributions of the flow variable S (77), velocity, 
and density of the stellar matter for a = 3 and various 
7e's are shown in Figures [1121 and [31 respectively. The 
same distributions but for a = 1.5 are shown in Figures 
141 5 1 and [6] Each line in these figures traces the time 
evolution of each fluid element after the arrival of the 
shock front (at 77 = 1). At the same time, each line in 
the range of < 77 < 1 represents the spatial distribution 
of the shocked fluid element before the shock breakout 
while lines with 77 < represents the spatial distribution 
after the shock breakout. The results show that the en- 
ergy variation at immediately behind the shock front pro- 
nounces the acceleration after the shock breakout rather 
than before the shock breakout (Figs. [H [5]). All the 
numerical results suggest that the velocities converge to 
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Fig. 4. — Same as Fig. [TJbut for a = 1.5 




Fig. 5. — Same as Fig. [2]but for a = 1.5 



constant values and that the density distributions con- 
verge to power-law profiles as t — * oo. Such asymptotic 
behaviors are discussed in next section. 

4.2. Asymptotic behavior 

From the numerical results, the derivative of the flow 
variable 5" is found to converge to a constant value, thus 
the asymptotic behavior of S should be 



S oc r\. 



(24) 



Using this relation, we derive the power-law exponents of 
velocity, pressure, and density of the ejecta with respect 
to the variable xq as 



a-2A+(l+A)7 
P OC Xq 



p«^ +A+1 . (25) 



The values of the exponents of p, p for various 7 C are 
shown in Tables [5] and [31 From these relations, we obtain 



p oc u 



-(a+A+l)/A 



(26) 



Thus the density distribution in the free expansion phase 
can be obtained by substituting u — x/t into this expres- 
sion. The exponents in this equation are also shown in 
Table |U 

The late time behavior for the case including the en- 
ergy variation is essentially the same as adiabatic case 
(7c = 7) except for the value of the exponents. 

4.3. Chapmann-Jouguet detonation 

We concentrate on the flow called "Chapmann- 
Jouguet detonation" that satisfies the Chapmann- 
Jouguet condition; the velocity of the gas relative to 
the shock front is equal to the local velocity of sound. 
iBarenblatt fc S ivashinskii (1970) showed that this con- 
dition can be expressed as the relation between 7 and 
7 e . We review the procedure and apply it to this study. 

From the Rankine-Hugoniot relations (JT])-©, the ve- 
locity of the gas relative to the shock front is 



\ui-U\ 



7e 



1 



7e + l 



-U, 



(27) 



Q. 

~3_ 




11 

Fig. 6. — Same as Fig. |3]but for a = 1.5 
the velocity of sound is 



V27(7c - 1) 




U. 



(28) 



7c + 1 

Equating the right hand sides of equations ([2T| and 
the condition is expressed as 

7c = 2 7 +l. (29) 

Using this expression, Tf and T s in (|A17[) and (|A20[) are 
found to take the same value. This means that the singu- 
lar point is lo cate d at the shock front. Then by equating 
R defined in (|A8|) at the shock front to that at the sin- 
gular point, we obtain the analytic expression for A as 

A=|. (30) 

On the other hand, the following relation is derived from 
equations (O and (fT5)l ; 



U oc p 



-A/a 



(31) 
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TABLE 2 

POWER-LAW EXPONENT OF p FOR 7 = 4/3, 



7e 

a t 11/3 3.0 2.0 1.5 4/3 1.2 1.1 1.01 1.001 1.0001 
3 1.000 1.095 1.417 1.719 1.886 2.078 2.303 2.758 2.925 2.978 
+00 3.667 3.698 3.806 3.906 3.962 4.026 4.101 4.253 4.308 4.326 
1.5 0.5000 0.5550 0.7000 0.8459 0.9275 1.023 1.136 1.370 1.459 1.488 
+00 2.500 2.518 2.567 2.615 2.642 2.674 2.712 2.790 2.820 2.829 



TABLE 3 

POWER-LAW EXPONENT OF p FOR 7 = 4/3, 



7e 

a t 11/3 3.0 2.0 1.5 4/3 1.2 1.1 1.01 1.001 1.0001 
3 3.000 3.000 3.000 3.000 3.000 3.000 3.000 3.000 3.000 3.000 
+00 5.000 4.953 4.792 4.641 4.557 4.461 4.348 4.121 4.037 4.011 
1.5 1.500 1.500 1.500 1.500 1.500 1.500 1.500 1.500 1.500 1.500 
+00 3.000 2.973 2.900 2.827 2.786 2.739 2.682 2.565 2.521 2.506 



In C-J detonation, substitution of the relation (|30|) yields 



U oc p 



1/3 



(32) 



which is independent of a. Therefore this relation can 
be applied to an arbitrary power-law density profile. For 
a shock wave without energy variation (7c = 7), substi- 
tution of the values in Table [T] yields 

which is also independent of a. 



(33) 



4.4. Energy loss limit 

Here we consider the situation that 7 e — ► 1. However 
the density diverges for the case that 7 C is exactly unity. 
Therefore introducing the parameter 6 as 

S = 7c - 1, (34) 

we treat the behavior of the eigen value A when 5 is 
close to zero. Fig. [7] shows the behavior in the range 
of 1(T 4 < 5 < 1CT 2 for 7 = 4/3, a = 3.0,1.5. In this 
range, the following empirical relations are derived by 
interpolation; 



A 



1.35 51 for a 
0.65(5 - 50 for a 



3.0, 
1.5. 



(35) 



These formulae suggest that when the shock becomes ra- 
diative, the shock velocity tends to be independent of the 
initial density. Thus the shock is no longer accelerated 
due to radiative energy loss and propagates at a constant 
velocity. 

5. ENERGY SPECTRUM 

In this section, we derive the energy spectrum from 
the asymptotic behavior of physical quantities described 
above. The energy spectrum is defined as the mass whose 
energy is larger than e; 



From the relation 
e is found to be 



M(> e) = J pdx. (36) 
(|2"5"|) and e = u 2 /2, the dependence on 
M(> e) oc e -(«+i)/ 2 \ (37) 



0.2 
0.1 
0.05 - 

0.02 - 
0.01 - 



0.005 



- 1 1 1 — 1 — 1 — 1 1 1 



T 1 1 1 1 — I I I I 
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10" 



10"' 

5 



10"' 



Fig. 7. — The behavior of A as a function of S for 7 = 4/3, a = 
3.0, 1.5 

The values of the power-law exponent for various 7 are 
shown in Table OH Here 7 is assumed to be 4/3. The re- 
sults for 7c = 27+I give the upper limit for the exponent, 
because 7 C takes the maximum value for the Chapmann- 
Jouguet detonation and then A yields the maximum. The 
exponents increase as the energy is deposited at the shock 
front and converge to the upper limit in the Chapmann- 
Jouguet detonation. On the other hand, the energy loss 
decreases the exponents and the exponents diverge at 
7c = 1 related to the situation when all the internal en- 
ergy is lost at the shock front. This result implies that the 
energy distribution of cosmic rays for the case 7 e > 7 (en- 
ergy deposition) becomes harder than that for the case 
7c = 7- 

6. CONCLUSIONS AND DISCUSSION 

In this paper, we have derived a self-similar solution 
to describe the propagation of the shock wave whose en- 
ergy varies at the shock front. Using the solution, we 
describe the propagation of the shock wave in a plane- 
parallel medium having a power-law density profile. The 
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TABLE 4 

VELOCITY SPECTRUM FOR 7 = 4/3, 



> 

a t 11/3 3.0 2.0 1.5 4/3 1.2 1.1 1.01 1.001 1.0001 
3 -3.000 -3.150 -3.790 -4.683 -5.384 -6.511 -8.611 -24.76 -80.37 -274.8 

+00 -5.000 -5.199 -6.054 -7.244 -8.178 -9.681 -12.48 -34.01 -108.2 -367.4 
1.5 -3.000 -3.174 -3.750 -4.586 -5.240 -6.288 -8.234 -23.05 -73.16 -245.4 

+00 -6.000 -6.291 -7.250 -8.644 -9.733 -11.48 -14.72 -39.42 -122.9 -409.9 



TABLE 5 

ENERGY SPECTRUM FOR 7 = 4/3, 













7c 












a 


11/3 


3.0 


2.0 


1.5 


4/3 


1.2 


1.1 


1.01 


1.001 


1.0001 


3 


-2.000 


-2.100 


-2.527 


-3.122 


-3.589 


-4.340 


-5.740 


-16.50 


-53.58 


-183.2 


1.5 


-2.500 


-2.645 


-3.125 


-3.822 


-4.367 


-5.141 


-6.862 


-19.21 


-60.96 


-204.5 



motion of the ejecta after the shock breakout is derived. 
The asymptotic behavior of the ejecta is also derived and 
then the energy spectrum in the free expansion phase is 
deduced. Considering the case for 7 C = 27+ 1, we derive 
the power-law exponent of the energy spectrum analyti- 
cally, which gives the upper limit for it. 

The calculation o f the radioactive energy input model 
for type la SNe fsee lColgate fc McKedl969ft showed that 
the explosion of an n — 3 polytrope leads to the outer 
density profile 

p oc u , (38) 

which corresponds to the case 7 C ~ 1.6 in this work. 
Then, the total energy per unit mass at the shock front 
becomes 



E*.at = 



1 



1 (39) 



,7c - 1 7- 1/ Pi 3 Pi 

and the energy source per unit mass becomes 

3 pi 

which is ~ 30% of the total energy at the front. For the 
Chapmann-Jouguet detonation, the outer density profile 
is 



(40) 



p oc u , 



which corresponds to 



and 



Etot ~ Y~p~{'- 



21 pi 

8 pi 



(41) 



(42) 



(43) 



This implies that ~ 78% of the total energy at the front 
is deposited by nuclear reactions. 

On the other hand, the shock wave in a SN II suffers 
from the radiative cooling. [Chevalier (1982]) showed that 
the outer density profile of the ejecta of SNe II is steeper 
than that of SNe I, and adopted the relation 

p(xu- 12 . (44) 



We assume that the progenitor of the SN II is a red 
giant, whose stellar envelope is convective (a = 1.5). 
The relation (|44f corresponds to the case 7 C w 1.18 in 
this work, which deduces to 



and 



Etot ~ 

23 Pi 

q ~ —z — , 

9 p{ 



(45) 



(46) 



then ~ 30% of the total energy is lost at the shock front 
by the radiative cooling. 

Here we consider some applicable limits of this solu- 
tion. First, from the comparisons with some numeri- 
cal simulations, the treatment of the energy variation in 
this self-similar solution is found to be valid. However, 
this solution can not determine the coefficients of phys- 
ical variables because we cannot express ki in equation 
(|18[) in terms of physical parameters specifying the phe- 
nomenon. As a consequence, we can not obtain absolute 
values for the physical variables. Second, we consider the 
behavior of the energy source q near the surface, where 
the total energy per unit mass diverges because its den- 
sity becomes zero. Simultaneously, the energy source di- 
verges to ±00 for the energy deposition/loss. We assume 
that the energy deposition is caused by some nuclear re- 
actions. Then the energy generation rate must decrease 
near the surface, because the density decreases while the 
shock front accelerates toward the surface. Therefore, 
this model becomes unrealistic near the surface. On the 
other hand, for the energy loss, the situation is realistic. 
As the total energy diverges, the energy source also di- 
verges to +00. This situation is understood as the radia- 
tive diffusion deprives the stellar matter of the diverging 
energy. 
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APPENDIX 
DERIVATION OF DISTRIBUTIONS 
Formulation 

From the basic equations of gas dynamics (fTTj) and (fT2|) . the equation of motion is expressed as 



fx 



1 (?£ 
P \dxj t 



1 



dp 
Po \dx 



Using the flow variable S(r]), the second derivative of x yields 

(fx _ x a „ 
dt 2 



(Al) 



(A2) 



Substitution of the relation into the right hand of equation (|A1|) yields 



1 



dp 
Po \dx 



7c "1 
7c + 1 V7c + 1 

7c -1 
+1 V7c + 1 

a - 2A 



„Q-2A 



dx a [S(r,) - (1 + \)r,S>(r,)P 



2 -2A-1 



k 2 x 



(A3) 



[[s - (i + x) v st 

Therefore, the equation of motion is rewritten in terms of S(r]) as 



7 (1 + A)r; 



2 Qll 



(i + xys 



7e 



a - 2A 



XS' + (1 + X)r,S" } 
[5 - (1 + A)^r +1 J ' 

AS' + (1 + X)t]S" 



[S-(l + X)r)S>] 



'p 



7(1 + Ah 



7c + 1 V7c + 1 

and we obtain the following ordinary differential equation; 

2( 7c - I)? (a - 2X)[S - (1 + A)r?S'] - 7 A(1 + X)r)S' 



[S - (1 + A)77S" 



17+1 



5" = 



(1 + A) 2 2 7 ( 7o - iy - ( 7c + l)7+i [5 _ (1 + A)?y 5/]7+i ' 
The distribution of physical variables are obtained by solving this equation under the initial conditions; 

S(l) = l, S'(l) 



(A4) 



(A5) 



(A6) 



(1 + A)( 7c + 1) 

Determination of eigen values 

Equation (|A5|) can be integrated to r/ — > — oo only for the eigen value A, because the derivative S" diverges at 
a singular point where its denominator becomes zero. In other words, the eigen value A is determined so that its 
numerator and denominator of the right hand side of (|A5|) vanish at the singular point. Setting the denominator of 
the left hand side of (|A5|) to be zero, it is found that the value of ij at the singular point is determined by the relation 



n : 



( 7c + i)7+n 1/2 



27(7c - I) 7 



[S- (1 + A)^'] ( 7 +1 )/ 2 . 



(A7) 



Therefore r\ takes positive value at the singular point. We determine the eigen value A as follows. 
Introducing new variables R and T as 

r) y(i+i) T = S-(l + X)r ] S' 1 



(A8) 
(A9) 

equation (|A5[) is converted to two first order differential equations by the following steps. First S' is rewritten with R 
and T as 

1 



S' = 



Then its derivative is 



S" 



1 + A 



1 + A 



-27/(1+7) 



(l-7)/(l+7) 



(R-T). 



7 



1+7 



(R-T) 



dR 



dT 



Substitution of equations (|A9|) - (|A1 1|) into equation (|A5|) yields 



(1 + A) 



1-7 
1 + 7 



(R-T) 



dR 



dT 



dlnrj d\nrj 



= 2 (7e - 1) 



d In 77 d In 77 



[a + (7 - 2)A]T - 7 Ai? 
2 7 ( 7 e-l) 7 -[(7e + l)T]7+i' 



(A10) 
(All) 

(A12) 
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On the other hand, the derivative of equation (|A8|) is 



S' = r] 



(1-7)7(1+7) 



1+7 



-R 



dR 
dlnrj 



Combining this expression with equation (|A10[) . we obtain 



(1 + A) 



dR 
dlnrj 



2(1 + A) 
1 + 7 



and substitution of this equation into equation (|Al2j) yields 



(1 + A) 



dT 



-XR 



2+ (1 - 7 )A 



T - 2 ( 7e - 1)- 



R-T, 



[a + (7 - 2)A]T - 7 Ai? 



'rflnT? 1 + 7 * 2 7 (7 c -l) 7 -[( 7o + l)T]7+i' 

Dividing equation (|A15[) by equation (|A14[) . we obtain the following ordinary differential equation; 



1 



2 + (l-7)A [ a + ( 7 -2)A]r- 7 Ai? 

1 + 7 J+Zl7c ij 2 7( 7c -ir-[( 7o + i)r]T+i 

2(1 + A) 



(A13) 



(A14) 



(A15) 



(A16) 



1 + 7 



R-T 



The initial conditions for the variables R and T are derived from its definition (|A8|) 

7e-l 



i?f = 1, If 



7e 



1 



9|, and the condition (| A6|) as 
(A17) 



For an arbitrary A, the integration of equation (|A16[) from the values of the shock front (Rf , T{ ) can not necessarily reach 
the other boundary, because the derivative dT/dR diverges . To avoid this situation, we require that the denominator 
and the numerator of equation (|Al6j) vanish at the same point, which leads to the following conditions 



2l 



[a+( 7 -2)A]T- 7 Ai? = 0, 



-(l + 7c )T 7+1 = 0. 



7c + 1 



(A18) 
(A19) 



Therefore, we determine the eigen value A so that equation (|A16[) can be integrated from (Rf,Tf) to (i? s ,T s ) given by 



R* = 



a + (7 - 2)A 
7 A 



27 fie - 1 



7o + l V7c + 1 



71 1/(7+1) 



T s = 



27 /7c-l 



7c + 1 V7c + 1 



1/(7+1) 



(A20) 
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